Shadow band in the one— dimensional large U Hubbard model 
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We show that the factorized wave-function of Ogata and 
Shiba can be used to calculate the k dependent spectral func- 
tions of the one-dimensional, infinite U Hubbard model, and 
of some extensions to finite U. The resulting spectral function 
is remarkably rich: In addition to low energy features typical 
of Luttinger liquids, there is a well defined band, which we 
identify as the shadow band resulting from 2kF spin fluctua- 
tions. This band should be detectable experimentally because 
its intensity is comparable to that of the main band for a large 
range of momenta. 
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The calculation of the spectral functions of models 
of correlated electrons is one the most challenging and 
largely unsolved issues of condensed matter theory. Al- 
though a number of numerical techniques can be used, 
e.g. exact diagonalization of finite clusters Q or quantum 
Monte Carlo simulations exact results are available 
only in very special cases, mostly for one-dimensional 
spin models |^ . As far as one-dimensional electron mod- 
els are concerned, most of the well established results 
have been obtained in the framework of the Luttinger 
liquid theory 0,^,||j7j , which is believed to be the correct 
description of the low energy properties of a large class 
of Hamiltonians. However, an accurate determination of 
the dynamical properties for all frequencies is so far still 
lacking. 

In this paper we perform such a calculation for the 
following one-dimensional models: 

i) The Hubbard model defined by the Hamiltonian 



(1) 



in the infinite U limit, which is also equivalent to the 
J — > limit of the standard t — J model; 

ii) An extension of the t — J model first proposed by 
Xiang and d'Ambrumenil S defined by the Hamiltonian 



^ = -^E(4.2,m^-f /i.e. 
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spin between them. The motivation to study this model 
is that, unlike the infinite U Hubbard model, there is 
an energy J associated to spin fluctuations, and this will 
give us useful indications about the 1/U corrections in 
the case of the finite U Hubbard model. 

Although the Hamiltonians of the two models are dif- 
ferent, they share the remarkable property that in both 
cases the eigenstates can be factorized flplEolO] as 



|/,iV) = |V^^(Q,{/}))®|xM(Q,/Q)) 



(3) 



where |^^(Q,{/})) 



is an eigenfunction of N non- 
interacting spinless fermions on L sites with momenta 



27r/j + Q (Ij are integers, j = 1 . . . N) and 
Ixa'(Qi/q)) is an eigenfunction of the one dimensional 
spin-i Heisenberg model with TV spins (we choose N 
as even integer not multiple of four) and momentum 
Q — 2ttJ/N, J integer. This momentum imposes a 
twisted boundary condition with phase e*'^ to the spin- 
less fermions |12 . For more details, see Ref. jl^. This 
wave-function has already been used by Ogata and Shiba 
to calculate the momentum distribution function and 
by Penc, Mila and Shiba to calculate the local spectral 
function of the infinite U Hubbard model . 

In the following, we will determine the full momentum 
dependence of the photoemission and inverse photoemis- 
sion spectral functions defined by 

Aik,u;) = Y^\{f,N+l\cljO,N)\'s{u;-Ef+' + E^) , 

B{k, c^) = E I (/' ^-l|Cfc,.|0> N)\' S{u;~E^ + Ef-') . 

As a result of the factorized form of the wave functions, 
the spectral functions can be obtained as a convolution: 

Uj' ,Q,(7 

B{k,co)^ D,iQ,u;')BQ{k,co~Lo'). (4) 



W ,U,(7 



where c are the usual projected operators and Vij ~ 

Ylj~2i{^~^i+j') ill tii6 exchange part of the Hamiltonian AQ(fc, — ^E] 
ensures that two spins interact as long as there is no other {/} 



A similar expression holds for the spectral function in 
the upper Hubbard band A™^{k,uj « [/), which 
we will not discuss here. In these expressions, AQ{k,uj) 
and BQ{k,oj) involve only the spinless fcrmion part of 
the wave function and are defined as 
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BQ[k,u:) = |«-i(g, {/})|6ol<'r)f 
{1} 

xSiLu-E^+Ef~')S{k~Pf-'+P^) , (5) 

where the momentum and energy of the states are given 
by P^' = eJi/cj and E'^' = -2tJ2fliCOskj, and 
where b and are spinless fermion operators. C„{Q,uj) 
and Dc{Q, oj) depend on the spin wave fmiction only and 
are given by 

a(Q,^) = ^|(x^+i(g,/Q)l4Vlx?)f 

/q 

x5{u - Ef+^ + E^) , 

D„{Q,uj)=Y,\{XN-i{Q.fQ)\ZoAx%'')^ 
Iq 

xSiu-E^ +Ef-'), (6) 

where Zq ^ appends a spin a to the beginning of the spin 
wave function \xn) making it iV + 1 sites long, and u 
is the hermitian conjugate of Zq ^. 

To evaluate the charge contribution, one needs matrix 
elements between states with different boundary condi- 
tions (6**2 for the final state, for the ground state 
p3[ ) . For Q ^ TT the overall phase shift {Q — tt)/L due 
to momentum transfer Q — tt to the spin degrees of free- 
dom gives rise to Anderson's orthogonality catastrophe 
[ p^ and the matrix elements n +i{Q , {I})\b\\i^'^^)\^ 
can be shown to be equal to 

L-^^-icos^^^TTsin"^^ 

j>i 

k' - k'- y h- 

xnsin^^n--"'^' (7) 

where kj (kj) are wave vectors with phase shift Q/L 
(tt/L). The restriction imposed by 6{k - Pf^^ + P^) 
is then implemented by restricting the sum over {/} to 
states which have the correct momentum. 

The calculation of the spin contribution is based on 
the spin-i Heisenberg Hamiltonian with N' sites 

N' 

Hspin = X! X! {^i^i+l - \^a,z) , (8) 

1=1 a=x,y,z 

with N' — N for the ground-state and TV ± 1 for the final 
states. The model of Eq. (H) corresponds to J" = J". 
For the infinite U Hubbard model, one has to consider 
the isotropic case and to take the limit J ^ 0. In that 
case, there is no energy associated to spin excitations, 
and we can write Co- (Q, = Ca{Q)5{uj) and D„{Q,uj) = 



Dcr{Q)S{uj). The functions Ca{Q) and D^iQ) have al- 
ready been studied in our previous paper . They can 
be calculated numerically with exact digonalizations (up 
to 26 sites) or with DMRG (up to 130 sites, keeping 
300 states per block). It turns out that there is a very 
strong singularity {n/2 — Q)^^^^ for Q < n/2 and some 
background coming from the higher order excitation tow- 
ers for Q > tt/2 in the case of Da{Q,uj). For Ca{Q,tjj) 
the situation is reversed and both are symmetric with 
respect to Q = 0. 

Using these results, it is straightforward to get the 
spectral functions for the infinite U Hubbard model. One 
just has to generate the quantum numbers for the charge 
part, calculate the corresponding energy, momentum and 
matrix elements, and perform the convolution in Q. The 
results are presented in Fig. |l| and || for a quarter- 
filled system. There are several interesting features to 
notice. In the low energy region near kp we can identify 
three structures. For k < kp there are divergences at 
Lu — Uc{k — kp) and uj — and a lot of spectral weight 
between them (peaks 'b' and 'c' on Fig. ^ ). There is 
also a small weight ('e') appearing on the other side of 
the Fermi energy for lu > —Uc{k — kp). For k > kp the 
spectrum is symmetric with respect to kp. If we remem- 
ber that the spin velocity Us vanishes for the infinite U 
Hubbard model, all these features are consistent with the 
Luttinger liquid calculations of Meden and Schonhammer 
and of Voit [Q. The small peak 'g' comes from higher 
harmonics. The dispersion of the charge part ('b') is ex- 
actly given by E{k) = —2tcos{\k\ + kp), in agreement 
with the observation of Preuss et al Q based on Monte 
Carlo results for U/t — 4. 

However, the Luttinger liquid picture does not ex- 
haust the features of the spectral function of Fig. 1 
and 2. For larger energies, or away from kp, there is 
a well defined band-like structure ('a') with consider- 
able spectral weight and a dispersion given by E{k) ~ 
— 2/:cos(— |fc|-ffcj?). We interpret this feature as a shadow 
band ITtl coming from the spin fluctuations which di- 
verge at 2kp. The scattering of the charges by these 
fluctuations produces an image of the main spectrum at 
k 2k p. This is very similar to the mechanism of the 
shadow bands proposed for the two dimensional model 
with strong antiferromagnetic fluctuations. This shadow 
band is responsible for the singularity at Zkp present in 
the momentum distribution function ]9| Jl^ , |l^ . Finally, 
there is a Van Hove singularity at ±2t which gives rise 
to a clear peak for wave vectors close to the extrema of 
the bands ('f'). 

Let us now turn to the model of Eq. (^). To get the 
spectral function, we need C^iQ, uj) and D„{Q, uj) for the 
Heisenberg model. This can be done numerically for the 
isotropic case [J^'^'^ = J) using Lanczos diagonalization 
of smaU clusters or DMRG |2^. We find that Co(Q, w) 
is zero for oj < — Jln2 -|- Ua\ sin(fc — ^)|, where Ua = f J 
is the spin velocity in the squeezed system, that it has 
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an inverse square root singularity at Q = 7r/2, and that 
the largest contributions come from the lower edge of 
the excitation spectrum. The main difference with the 
infinite U case is that the spin fluctuations have an energy 
of order J, so that the spin velocity Us = u„L/N does not 
vanish anymore. The low energy part of the spectrum has 
then exactly the form predicted by the Luttinger liquid 
theory. 

For the XY case (J^'^ — J, J'^ = 0) one can give a 
closed expression for Ca{Q^uj) and D„[Q^uj) after map- 
ping the problem onto non-interacting spinless fermions 
by a Jordan- Wigner transformation. After some algebra, 
the matrix elements |(xAr+i(Q, /q)|4,^|x^^)P of Eq. (|) 
can be obtained as 

M / 

i=l ]>i 

xJIsin^^n--"^' (9) 

j>i id 

where qj and q'^ are the momenta of the M spinless 
fermions representing the —a spins on the N and iV + 1 
site lattice. They are quantized according to q' = 
2ttJj/{N + 1) and q = 2-Kjj/N, where J, and Jj are 

integer quantum numbers, and fq = {Jpj = l.-M}. 
The total momentum and energy of \xN+iiQi fq)) are 
given by Q = J2j1j = JJ2j^osq'y De- 

tails will be given elsewhere |Q. A similar expression 
holds in the case of D^{Q,uj). This formulation also al- 
lows one to derive analytical results. For instance, the 
static function a;(0 — > j,a) introduced by Ogata and 
Shiba |] can be shown to have the asymptotic behav- 
ior cx j~^/^cos(^j -I- j). Thanks to this mapping, one 
can calculate the spectral function with the same accu- 
racy as for the infinite U Hubbard model. The results 
are shown in Fig. ^ for a quarter-filled system. It is 
essentially the same as that of the Hubbard model, ex- 
cept that at low energies an extra peak 'd' accounting 
for the extra exponents in the spin part of the XY model 
has appeared (this peak has nothing to do with peak 'g' 
on Fig. 2.) Due to finite J, both 'c' and 'd' follow the 
Lu = Ucr COS § dispersion. Furthermore, we can see that 
the shadow band ('a') and the Van Hove like singularity 
('f ') are broadened by the spin fluctuations. 

Finally, let us comment on the experimental implica- 
tions of the present results. It would be most interesting 
to observe the shadow band in angular-resolved photoe- 
mission or inverse photocmission experiments on quasi- 
one dimensional conductors. The intensity of that band 
in the previous calculations is certainly big enough for it 
to be detected. What about the experimentally more rel- 
evant case of the Hubbard model with finite U? In that 
case the factorized wave functions are no longer eigen- 
functions of the Hubbard model, and there are two types 
of 1/C/ corrections to the spectral functions. The first 



type is due to the energy coming from the spin part with 
an effective coupling J ^ ^{n — ^"^^J" ) ■ expect 
these corrections to be very similar to those of the model 
of Eq. (H), and the main effect is to give a finite veloc- 
ity to the spin excitations. However, there are also 1/U 
corrections entering the matrix elements of the spinless 
part of the wave-function. We can anticipate that they 
will have two effects on the spectral function. They will 
produce a transfer of spectral weight to the upper Hub- 
bard band which, according to Eskes and Oles |2^, will 
be small except very close to half-filling, and they will 
modify the power laws of the singularities. So, at least 
not too close to half-filling, the shadow band seems to 
be robust against 1/U corrections. Whether this remains 
true for small values of U is not clear yet. Let us just 
mention that, according to recent numerical results ob- 
tained by Maekawa et al. ||2^ in a study of the spectral 
function of the Hubbard model for U/t — 10 based on 
Lanczos diagonalization of finite clusters, there seems to 
be a structure in addition to the Luttinger liquid fea- 
tures, suggesting that U/t ~ 10 is already large enough 
to guarantee the presence of a well defined shadow band. 

We thank to Y. Kuramoto, H. Fukuyama, M. Imada, 
D. Poilblanc, K. Vladar and A. Zawadowski for useful 
discussions. 
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FIG. 1. One particle spectral functions of the U — > +00 
Hubbard model for L = 228 sites and A'^ = 114 electrons with 
Fermi momentum kp = 7r/4. 

FIG. 2. The same as Fig. 1, but for some selected momenta. 
Some parts of the spectra are multiplied by 10 and are shown 
with dashed lines. 

FIG. 3. Spectral function for the model of Xiang and 
d'Ambrumenil with XY exchange, J = 0.4i, L = 228, 
N = 114 and ef = — J/tt. Some parts of the spectra are 
multiplied by 10 and are shown with dashed lines. 
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